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1-^ ■ We present new interaction potential curves, calculated from first-principle, for the 

He(^S)- ■ ■ H2 and He(^S)- ■ ■ Ar systems, relevant in Penning's ionization experiments. 

^ • Two different approaches were applied: supermolecular using coupled cluster the- 

. ! ory (CC) and perturbational within symmetry-adapted perturbation theory (SAPT). 

.^ ' Both methods gave consistent results and the potentials were used to determine the 

c/3 ! 

^' positions of shape resonances in low collision energy scattering regime. We found a 

^ ■ 
, P^ i good agreement with the most recent scattering experiment of Henson et al. [Science 

338, 234 (2012)]. 

> ■ 

^ . In addition, we investigated two other dimers, composed of metastable Ne and 

00 ' 

\Q ' ground state He and Ar atoms. For the Ne( P)---He system a good agreement 

.^ ■ between CC and SAPT approaches was obtained. The Ne( P)- ■ ■ Ar dinier was de- 

Z^ ', scribed only with SAPT, as CC gave divergent results. Ne* systems exhibit extremely 

small electronic orbital angular momentum anisotropy of the potentials. We attribute 

kJ • this effect to the screening of the open 2p shell by the singly occupied 3s shell. 
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I. INTRODUCTION 

The quest for taking control over chemical reactions and scattering processes by manipu- 
lating the movement of molecules led in the last decade to a rapid developement of methods 
which use the electrioi"- or magnetic fields^>^ control combined with molecular beam tech- 
niques.- The possibility for quantum control opens in the low-energy regime, when the 
kinetic energy of the molecule becomes comparable to the perturbations due to the external 
electromagnetic field.- Under these conditions the scattering is dominated by several partial 
waves only, facilitating studies of such phenomena as scattering resonances. Moreover, elab- 
orate detection techniques allow to study the state-resolved cross sections and the propensity 
rules, governing the distribution of the inelastic scattering products.- Finally, the experi- 
ments on low-energy scattering provide unique opportunities for verification of the potential 
energy surfaces of colliding species, which is particularly beneficial for developing theory of 
intermolecular forces. 

Recent progress in this field involved theoretical works in the scattering theory^ and ex- 
periments employing metastable noble gas atoms interacting with other species.— »ii Studies 
of Henson et alM- revived the interest in Penning ionization reaction (PI) - in their work a 
beam of metastable helium atoms in the 2'^S state, with their velocity controlled by a time- 
varying magnetic field, has been merged with a supersonic beam seeded with H2 molecules 
or Ar atoms. Thanks to the velocity control, for the first time it was possible to observe 
sub-Kelvin scattering resonances with sub-milikelvin resolution. 

Understanding of PI is very important since this process occurs in plasma^^ii^ and surface 
chemistry.— The reaction takes place when molecules collide with species excited to a very 
high energies, such as metastable 2^S He or ^P Ne atoms. Provided that an atom is excited 
above the ionization threshold of the molecule, the collision quenches the atom down to the 
ground state releasing internal energy, which in turn can ionize the molecule: A +B — y A + 
B+ + e-. 

Metastable He and Ne atoms have also been of great interest to cold matter physics, 
in particular the Bose-Einstein condensate of metastable He atoms has been obtained.— li^ 
It was demonstrated that light atoms at very low temperatures could be used for efficient 
sympathetic cooling of moleculea^'i^ and metastable helium atoms in fiK regime could be 
considered as coolants, provided that the anisotropy of the potential energy surface is not 



very large.— i^ 

Because of the lack of interest, the theoretical effort to improve and develop potentials 
for metastable atoms interacting with molecules has been abandoned for more than 10 
years. However, new perspectives offered by the experiments with metastable atoms require 
knowledge of accurate potential energy surfaces of those systems. 

There is a principal difficulty in treating electronic states of the A*-B type that are 
submerged in the continuum of states of the (AB)"'"+e type: the variational principle, which 
is the fundamental tenet of quantum chemistry, is expected to drive such states down either 
to the ground state, or one of the excited states or - because of the coupling with continuum 
states - to some delocalized state corresponding to the fragmentation into an ion, a molecule, 
and a free electron. Thus, during the wave function optimization one should seek for a 
constraint minimum of the energy functional that does not obey the Aufbau rule. 

Interaction potentials for systems undergoing PI have been obtained with either ab initio 
or semi-empirical methods. Fully ab initio approaches encompass the "stabilization method" 



(e.g., see Ref. |20[ ). the configuration interaction approach with embedded Feshbach projector 
operators^! and the Siegert-eigenvalues techniques.— It is, nevertheless, the semi-empirical 
method based on fitting parameters of model potentials to the experimental results that 
enabled construction of the most accurate PI interaction potentials. Different models have 
been suitably adapted over the years: a one-electron model potential by Siskai^, Lennard- 
Jones potentials and so called piecewise potentials, constructed with components insuring 
a proper short- and long-rage behaviour. Ultimately, in order to describe He*(2 ^S)- ■ -112 
interaction Martin and Siska introduced a non-piecewise potential surface function based on 
a modified Tang-Toennies model2^ (for a thorough review of previous theoretical works see 



Ref. 



251). 



The most straightforward among many-electron ab initio methods is the "stabilization 
method" (SM). The wave function is expanded in an atomic/molecular basis set, followed 
by diagonalization of the Hamiltonian matrix, thus all standard procedures designed for 
electronic bound states may be applied. An obvious disadvantage of SM is a possibility of 
obtaining continuum-like solutions in close proximity to the true resonances. Fortunately, 
such spurious solutions may often be identified due to their strong basis-set dependence. 

In this work we present state-of-the art calculations of intermolecular potentials for He*- 
and Ne*-molecule dimers. To this end, we applied two independent methods: (i) the super- 



molecular approach based on the coupled cluster method with judiciously prepared reference 
states, and (ii) the perturbation approach involving symmetry-adapted perturbation theory 
(SAPT) for open-shell species.— 

Previous works exploring application of different formulations and levels of SAPT ex- 
pansion to excited-state systems are sparse: they can be found in Refs. l27H3ll. Interactions 
of a metastable helium atom were addressed in two of them - in Ref. |30| Przybytek et al. 
performed a detailed study of the convergence behaviour of SAPT expansion for He('^S)- ■ ■ H 
interaction, whereas Zuchowski et alM- calculated the dispersion energy contribution in the 
random-phase approximation for the He(^S) dinier. 

In the case of PI reaction, the SAPT methodology is appealing, as it ensures a fixed 
number of electrons at each monomer. The issue of concern, however, is the convergence of 
SAPT, as shown e.g. for dimers of metals.— 12^ 

The systems which we investigate in this paper either undergo Penning ionization: 
He*---H2, He* ■ ■ ■ Ar (studied also by Henson et al.—), and Ne*---Ar— , or provide an 
auxiliary model for testing SAPT performance for such systems, as the Ne* • • • He dinier. 

II. COMPUTATIONAL DETAILS 

For both complexes of metastable helium we were able to perform supermolecular unre- 
stricted coupled cluster singles and doubles with noniterative triples correction (UCCSD(T)) 
calculations.— It should be stressed, that convergence of CCSD(T) for excited states calcula- 
tions cannot be taken for granted. Therefore, we performed them carefully by the following 
procedure. First, we ran separate monomer calculations where desired states of monomers 
were set: in our case, one monomer was a metastable helium, the other was either the 
ground state H2 molecule or an Ar atom. The resulting orbitals were then used to build 
the starting wave function of the dinier as a product of two monomer functions. During 
the optimization of the ROHF wave function we did not use of the Aufbau principle, and 
in each iteration of ROHF, orbitals were reordered to obtain maximum overlap with closed- 
shell/open-shell/virtual spaces from the previous iteration. For the Ne* ■ ■ ■ He system the 
UCCSD(T) calculations were straightforward since the state is the first triplet state. The 
UCCSD(T) results subsequently served as a benchmark for the performance of SAPT cal- 
culations. It should be noted, that the CC approach did fail for the potential curves of the 



Ne* ■ ■ • Ar system. 

Both metastable He('^S) and Ne(^P) are high-spin open-shell atoms and can be treated 
with SAPT based on restricted open-shell or unrestricted Hartree-Fock and Kohn-Sham 
description of the monomers, SAPT(ROHF) and SAPT(ROKS}^ or SAPT(UHF) and 
SAPT(UKS)^, respectively. The method is based on the expansion of the interaction energy 
in the power series of the intermolecular interaction operator, which up to the second order 
takes the form: 



pSAPT _ p(l) , p(l) , p(2) p(2) p(2) p(2) /-,^ 

-^int — -^elst "T -^cxch "•" -^ind "•" -^disp "•" -^cxch-ind "•" -^exch-disp' V / 



where on the right-hand side the consecutive terms represent the electrostatic, first-order 
exchange energy, induction and dispersion energies, and second-order exchange- induction 
and exchange-dispersion energies. In this work we choose the unrestricted, SAPT(UHF) 
and SAPT(UKS) formalisms. The second-order terms have been calculated at the coupled 
level of theory.— In SAPT(UKS) calculations the PBEO xc functional^^i^i was asymptotically 
corrected^^, with vertical ionization potentials calculated using a difference method. 

For UCCSD(T) calculations double augmented, correlation consistent quadrupole-zeta 
basis set (d-aug-cc-pVQZ) of Dunning et ai— was chosen for helium and hydrogen atoms. 
For Ar we used aug-cc-pVQZ basis set. Ne('^P)---He complex was described with d-aug- 
cc-pVQZ basis supplemented with 3s3p2d2f bond functions (with exponents equal to: sp 
0.9, 0.3, 0.1, df 0.6, 0.2) placed at the midpoint of the intermolecular distance. In SAPT 
calculations different basis sets had to be chosen in order to insure convergence of the 
metastable monomers. For a metastable helium, d3Z basis set of Przybytek^ was selected. 
Hydrogen was described by the aug-cc-pVQZ basis. In the case of He(^S)- ■ ■ Ar complex we 
chose the aug-cc-pVQZ basis set for argon and a set of 3s3p2d bond functions was added 
in order to saturate the dispersion interaction. In Ne('^P)---Ar aug-cc-pV5Z basis was 
extended with a set of additional diffuse functions of each angular symmetry, according to 
the even-tempered scheme implemented in MOLPRO.— 
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FIG. 1. a) Isotropic (Vq) and leading anisotropic terms (V2) of the He(^S)---H2 interaction, b) 
interaction potential for He(^S)---Ar calculated at different levels of theory. See the text for 
description of the calculation details. 

III. NUMERICAL RESULTS 

A. He(^S)- ■ ■ X systems 



Fig. [T] shows calculated potential energy curves for the He*-X systems and the param- 
eters of the interaction potentials, Re, De and Cq, are given in Table HI In the case of 
the He(^S)- ■ -112 system (Fig. [T|l) we have shown the isotropic (Vq) part of the interaction 
potential along with the leading anisotropic terms (^2)- 

We have investigated both SAPT(UHF) as well as SAPT(UKS) approaches as described 



in Ref. 



26 



5|. Analysis of the Vq potential parameters for the He(^S)---H2 system reveals 
that both SAPT schemes remain in a reasonably good agreement with UCCSD(T) reference 
values. The good performance of SAPT(UHF) can be attributed to the basis set optimization 



procedure used in Ref. |40|. Moreover, for a two-electron helium atom in its ^S state we 
found SAPT(UHF) to give a better description of electron spin density than SAPT(UKS). 
The erroneous behaviour of the PBEO spin density is reflected, for instance, by the static 
polarizability DFT calculations lead to the value of 258.79, whereas UHF predicts 317.04 
a.u, much closer to the benchmark result of 315.63 a.u..— 

UCCSD(T) places the minimum at 10.8 Oq with the depth of 9.65 cm~^, whereas 
SAPT(UHF) leads to well depth of 7.88 cm~^ at larger distance of 11.4 oq, which amounts 



6 



to 18% difference. These results are in good agreement with the potential derived from 
the experiment of Henson et al.—, where the well depth of 8 cm~^ has been found at R^ 
of approximately 11.5 Uq. For comparison, the most recent semi-empirical potential from 



Refs. 



24 and 



43l predicted the minimum at a similar distance, however its well depth twice 
as shallow (See Table H]). 

The leading anisotropy of the potential energy surface, V2 is very small. In the van der 
Waals minimum region both methods agree very well, however, the UCCSD(T) potential 
anisotropy is much larger in the short-range repulsive region. 

It is also interesting to compare results for the global minimum of the potential energy 
surface. Both methods give the minima at the linear geometry: UCCSD(T) predicts a well 
depth of 12.86 cm~^ at 10.6 a^, while SAPT leads to shallower minima. SAPT(UHF) gives 
a well depth of 9.28 cm^^ at 11.4 uq and SAPT(UKS) gives minimum of 9.23 cm^^ at 11.2 
ao- 

Similar to the He* ■ ■ ■ H2 system, the well depth of the He* ■ ■ ■ Ar potential energy curve 
is consistent for all of the applied methods (within 7%). The equilibrium distance of the 
SAPT potential curve is slightly larger than in case of CCSD(T) which is consistent with 
the behaviour of Vq of the He(^S)- ■ ■ H2 potential. 

In Tab. [J we present the (isotropic) Cq coefficients. Apart from comparing the calculated 
Cg coefficients with reference values found in the literature, we have applied the Tang's 
formula^: 

(2) 



a 
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6.AS 
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b) *-^6,Ayl 



where a^ stands for static polarizabilities and Cq^xy denote Cq coefficients for homo- and 
heteronuclear dimers. Static polarizabilities for He* (a^g, = 315.63 a.u.) and H2 (oi^^ = 

respectively. For argon atom we used a^^ = 



5.18 a.u.) were taken from Refs. 



42 



and 



45 



11.078 a.u. and C^^ArAr = 64.42 a.u. values from Ref. |46|. Finally, C^^He'Re* = 3276.67 a.u. 
and 6*6,^2^2 = 12.058 a.u. were found in Refs. |47| and ? , respectively. As summarized in 
Table [T] all of the applied methods lead to consistent results for Cr coefficients and stay in 



48 



and 



49|) 



agreement with previous theoretical prediction (see Refs. 

In order to test the He* ■ ■ ■ H2 and He* ■ ■ ■ Ar potentials we calculated the position of 
the shape resonances in the elastic low-energy (up to 1 meV) scattering of He(^S) with the 
molecular hydrogen and argon. Since the ah initio methods applied in this paper provide only 



the description of the entrance channel of PI process, we are not in a position to predict the 
reaction rates. Nevertheless, the enhancement of the PI reaction rates attributed by Henson 
et ai— to shape resonances is entirely the entrance- channel effect, thus they correspond 
exactly to the calculated elastic scattering cross sections in the reactive scattering. 



Prediction of resonances is a particularly demanding test of the quality of the potential 
curve. The calculated elastic cross sections for He* ■ ■ ■ H2 and He* ■ ■ ■ Ar systems are shown 
in the Fig. [2] and the overall pattern of the resonances is nearly identical for UCCSD(T) and 
SAPT potentials. For the He* ■ ■ ■ H2 system, for which we used only Vq potential in scattering 
calculations, we found strong shape-resonances: the first near 0.08 meV for UCCSD(T) and 
0.1 nieV for SAPT, and second at 0.3 meV for both potentials. They originate from the 
partial wave components with end-over-end quantum number I = 3 and 4, respectively. 
In the experiment both resonances were found at lower energies, ca. 0.02 meV and 0.2 
meV, respectively. Note that Henson et alM- attributed those resonances to / = 5 and 6, 
respectively, whereas in our work they appear for / = 3 and 4. In contrast to what is shown 
by Henson et ai, our adiabatic curve associated with the / = 6 partial wave exhibits no 
attractive minimum. 



In the case of He* ■ ■ ■ Ar system, both UCCSD(T) and SAPT potentials lead to slight shifts 
in resonance positions. The position of the first shape resonance, at ca. 0.02 meV, is in a 
very good agreement with experimental findings.— The second and third strong resonances 
in the experimental studies were found at approximately 0.1 meV and 0.3 meV, respectively. 
These features are, again, properly reproduced by calculations with the potentials obtained 
in this paper (peaks at ca. 0.2 nieV and 0.4 nieV). 

8 



c 
o 

u 

O) 
l/l 

I/) 

o 
u 
u 

in 
— 



SAPT 
UCCSD(T) 




Collision Energy / meV 



FIG. 2. Shape resonances in the elastic cross sections for the scattering of He* • • • Ar (upper panel) 
and He* • • • H2 (lower panel) systems. The arrows depict the strongest resonances which we compare 
to the experimenliii. 



TABLE I. Parameters of the isotropic part of the interaction potential, Vq, for 
He('^S)- • • H2 and He('^S)- ■ ■ Ar interaction potential together with Cg coefficients. TF 
denotes results obtained with Tang's formula 



He(3S)---H2 



He(3S)---Ar 



Method 



Re [flo] De [cm ^] Ceo [a.u.] Re [ao] De [cm ^] Ce [a.u.] 



CCSD(T) 



10.80 9.65 97.^ 



10.09 27.31 214.1 



SAPT(UHF) 11.38 7.88 112.8 



10.75 25.68 220.0 



SAPT(UKS) 11.24 7.62 104.9 



10.61 25.63 206. 



TF. 



100.2 



216.4 



Ref. 



11.72^) 4.72^) 109.6"^) 



9.77") 41.27") 220.9^^) 



^) See Ref. M and M ^^ See Ref. M ") See Ref. M '^^ See Ref. |49 



B. Ne(3p)- ■ -X systems 

For the Ne('^P)- ■ ■ Ar complex we were not able to obtain convergent CCSD(T) results and 
SAPT seems to be the only method allowing for an insight into the nature of Ne(^P)- ■ ■ Ar 
interaction. In order to test the performance of SAPT in this case, we decided to choose a 
similar Ne(^P)- ■ ■ He system, in which PI does not occur. For both of the studied complexes 
the SAPT(UKS) results proved to be of practically the same accuracy as the SAPT(UHF) 
ones, hence, only SAPT(UHF) potentials are discussed. 

In Fig. |3^ and Table HT] we compare UCCSD(T) and SAPT(UHF) interaction potentials 
for ^S and '^11 states of Ne('^P)- ■ -He. The interaction potentials for that system are ex- 
tremely weak: the well depths are very shallow, on the order of 1 cm~^, and the minima 
occur at a very long range (> 12 ao). The performance of SAPT(UHF) in comparison with 
UCCSD(T) is similar to the He*- ■ ■ H2 case. The values of the potentials are very small 
due to strong cancellation of exchange and dispersion energies - thus, they are strongly 
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TABLE II. Global minima for Ne('^P)- • • He/Ar predicted at different levels of theory. 

Ne---He(2s) Ne- • • Ar (^S) 

Method Re [ao] De [cm"^] Ce [a.u.] Re [oq] De [cm^^] Cq [a.u.] 



CCSD(T) 


12.22 1.414 24.6 


— — — 


SAPT(UHF) 


13.25 0.998 29.1 


10.8 24.68 231.24 




Ne- • • He (^H) 


Ne- • • Ar (^n) 


CCSD(T) 


12.15 1.393 24.3 


— — — 


SAPT(UHF) 


13.14 0.980 28.3 


10.6 24.40 225.11 



susceptible to relative errors. SAPT(UHF) correctly recovers the S — 11 anisotropy in this 
system, which is around 0.03 cni~^ in the niininiuni and therefore should be regarded as a 
subtle effect. 

Fig. |3b and Table HTl present results for Ne(^P)- ■ • Ar. SAPT predicts a minimum at 10.6 
bohr for ^S and 10.8 bohr for "^11 state, 24.68 cm~^ and 24.40 cm~^ deep, respectively. Such 
extremely weak anisotropy of Ne* is, in fact, not suprising. A similar effect was observed 
for lanthanides and some transition metals in their ground state where it was attributed to 



51 



-53 



In the 



the effective screening of the doubly- filled outermost s shell, see, e.g., Refs. 
metastable neon atom such screening is even more pronounced. A weak anisotropy of Ne('^P) 
can be visualized in the following orbital picture: the electron lying in the outer 3s orbital 
effectively screens the inner 2p shell, which has much lower (r) value, causing substantial 
suppression of the anisotropy. It also manifests itself in a small quadrupole moment of the 
metastable Ne atom (0.092 a.u. calculated as expectation value with the CISD method), 
an order of magnitude smaller than quadrupole moments of the other second-row P-state 
atoms (e.g. 0.92 for O '^P atom). In their work Gregor and Siskai^^i also postulated that effect 
and therefore assumed a single triply degenerate S-II component in their model potential 
for Ne('^P)- ■ ■ Ar. They placed the minimum at 9.45 Qq with the well depth of 43.72 cm~^. 
The Cq = 225.11 a.u. obtained here with SAPT(UHF) agrees fairly well with the value of 
204.68 a.u. calculated using the theory of Proctor and Stwalley^ and given by Gregor and 
Siska.^ 
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FIG. 3. Potential energy curves for a) Ne('^P)- • • He, b) Ne(^P)- • • Ar calculated at different levels 
of theory. See the text for description of basis sets. 

IV. CONCLUSIONS 

We have investigated the potential energy curves of the He* ■ ■ ■ H2 (Vq and V2) and 
He* ■ ■ ■ Ar systems with two different methods: UCCSD(T) and SAPT(UHF/UKS). In both 
cases we have obtained a consistent picture of the interaction with minor differences in pre- 
dicted Dg and R^. We have used these potentials to calculate the positions of the shape 
resonances in low-energy scattering obtaining a good agreement with recently reported state- 
of-the-art scattering experiment of Henson et al..— 

We have also investigated two systems containing metastable Ne atoms: in case of 
Ne* ■ ■ ■ He system we have obtained a good agreement between UCCSD(T) and SAPT which 
suggests that the SAPT description of Ne*Ar dinier should be correct. We have revealed 
strikingly small electronic orbital angular momentum anisotropy of these potentials which 
can be explained by a strong screening of the open 2p shell by the outermost, singly occupied 
3s shell. 

As far as the performance of SAPT is concerned, we have exploited the unique feature 
of this method: fixing the number of electrons at each monomer to prevent an electron- 
hopping. This approach is expected to be general, provided that proper reference states of 
the monomers can be obtained (either Hartree-Fock or Kohn-Sham wave functions). An- 
other prerequisite of applying SAPT must be its convergent behaviour in the lowest orders. 
It cannot be taken for granted in general, for instance when both monomers are largely po- 
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larizable - in this last case, one could apply appropriately modified Pauli-Blockade method^ 
combined with second-order dispersion contribution obtained with SAPT. 

In the near future we hope to take advantage of SAPT to survey the potential energy 
surfaces of He*-paramagnetic molecule (or atom) systems exploring the possibility of their 
sympathetic cooling. 
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